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REAL-HELIUM HYPERSONIC FLOW PARAMETERS FOR PRESSURES 
TO 3600 ATMOSPHERES AND TEMPERATURES TO 15 000° K 

By Charles G. Miller HI and Sue E. Wilder 
Langley Research Center 

SUMMARY 

The virial form of the equation of state is used to derive a method for determining 
reservoir, nozzle -throat, free-stream, and normal- shock conditions for real-helium flow. 
This method is derived primarily to meet the needs of hotshot- type wind tunnels where 
the basic input data are measured reservoir pressure and stagnation pressure behind a 
normal shock. However, the method is derived to include an additional measurement of 
free-stream density, free-stream velocity, or reservoir temperature. The method is 
applied for reservoir pressures to 3600 atmospheres, temperatures to 15 000° K, and a 
range of free-stream Mach numbers from 10 to 100, and the results are presented in the 
form of nondimensionalized reservoir thermodynamic and nozzle-throat parameters and 
as simple correction factors which may be applied to ideal-helium flow parameters. 

These calculations are restricted to reservoir densities less than the critical density of 
helium and to theoretically uncondensed free-stream flow. 

INTRODUCTION 

Recent tests in the Langley hotshot tunnel have demonstrated the feasibility of 
operating this facility with helium as the test gas at Mach numbers from approximately 
30 to 80. These preliminary results were obtained for reservoir pressures from approxi- 
mately 200 to 1000 atmospheres and reservoir temperatures from 800° to 4000° K. For 
these reservoir conditions, appreciable errors can result from neglecting real- gas effects 
in determining the flow parameters (refs. 1 and 2). Therefore, a method employing the 
virial form of the equation of state was derived for determining reservoir, nozzle- throat, 
test- section free-stream, and test- section normal- shock conditions in real helium and is 
applicable for reservoir temperatures from 200° to 15 000° K. 

This method was derived to be used primarily with measured reservoir and test- 
section stagnation pressures as inputs and was based on the assumption of uniform density 
throughout the reservoir during a tunnel test. In the event that another measurement such 
as free-stream density, free-stream velocity, or reservoir temperature is made to avoid 



the assumption of uniform density throughout the reservoir, the method was also derived 
with the capability of utilizing an additional measurement as input. 

Real-helium correction factors (ratio of real-helium value of a particular flow 
parameter to its ideal-helium value) were calculated and presented in the form of working 
charts for reservoir pressures to 3600 atmospheres, reservoir temperatures to 15 000° K, 
and a range of free- stream Mach numbers from 10 to 100. A similar technique was used 
in reference 1 for which approximate results were obtained above 590° K using the Beattie- 
Bridgeman equation of state with constants which were determined from experimental data 
obtained at temperatures up to only 590° K. The present method employing the virial 
equation of state, for which the virial coefficients of helium were obtained from refer- 
ences 3 to 7 for temperatures up to 15 000° K, was used to obtain more accurate real- 
helium effects for temperatures above 590° K. The present calculations are restricted 
to reservoir densities less than the critical density of helium and to theoretically uncon- 
densed free- stream flow. 


SYMBOLS 


a speed of sound 

a r reference speed of sound 

A area 

B second virial coefficient 

C third virial coefficient 

Cp specific heat at constant pressure 

c v specific heat at constant volume 

d diameter 

D fourth virial coefficient 

h specific enthalpy 

M Mach number 
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p 

pressure 

q 

dynamic pressure 

R 

universal gas constant 

s 

specific entropy 

T 

temperature 

U 

velocity 

V 

reservoir volume 

Z 

compressibility factor 

r 

ratio of specific heats 

p 

coefficient of viscosity 

p 

density 

Per 

critical density 

T 

time 

Subscripts: 


c 

initial reservoir conditions 

calc 

calculated 

con 

condensed 

i 

ideal- gas conditions 


m 


measured 



reservoir stagnation conditions 


t,l 

t,2 stagnation conditions behind normal shock 

oo free- stream conditions 

2 static conditions immediately downstream of shock 

Superscripts: 

~ approximate value 

* nozzle-throat conditions 

BASIC THERMODYNAMIC RELATIONS 

The virial form of the equation of state is used in conjunction with various thermo- 
dynamic relations to determine expressions for enthalpy, entropy, speed of sound, and 
specific heats, each as a function of density and temperature. 

The virial form of the equation of state of a pure gas can be written as 

p = pRT[l + pB(T) + p 2 C(T) + p 3 D(T) + . . ] = pRTZ (1) 

where B(T), C(T), D(T), . . ., are the virial coefficients and are functions only of the 
temperature. The second virial coefficient B(T) represents the deviations from the 
ideal-gas state arising from interactions involving two molecules of the gas; the third 
virial coefficient C(T) represents deviation arising from interactions involving three 
molecules; and so forth. In the present study, interactions involving five or more atoms 
are neglected; hence, the compressibility factor Z was determined from 

Z = 1 + pB(T) + p 2 C(T) + p 3 D(T) (2) 

The second, third, and fourth virial coefficients of helium were determined by Amdur 
and Mason (ref. 7) for a temperature range of 1000° to 15 000° K. The interatomic force, 
laws used to determine the virial coefficients were derived from molecular-beam- 
scattering measurements. The virial coefficients of reference 7 for helium are 

B(T) ~ (1.3436 X 10" 2 )(15.8922 - In T) 3 + . . . (3) 

C(T) ~ (9.0263 X 10“ 5 )(15.8922 - In T) 6 (4) 
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and 


D(T) ~ (7.0341 X 10-7)(15.8922 - In T) 9 (5) 

where B(T), C(T), and D(T) are in cm 3 / mole, (cm 3 / mole) 2 , and (cm 3 /mole) 3 , respec- 
tively, for T in °K. The contribution to B(T) of the second and subsequent terms of 
equation (3) as presented in reference 7 is negligible and is therefore omitted. As shown 
in figure 1, where the second virial coefficient is plotted as a function of temperature, the 
values of B(T) obtained from equation (3) for temperatures less than about 1000° K dif- 
fer from the experimental results of references 3 to 6, with this discrepancy decreasing 
with increasing temperature. The second virial coefficient is negative at very low tem- 
peratures as a result of the atomic attraction forces, becomes positive at somewhat higher 
temperatures as repulsion forces start to dominate, passes through a maximum, and then 
decreases with further increase in temperature as the repulsion forces dominate. There- 
fore, this discrepancy between the calculated values of B(T) from equation (3) as 
obtained from reference 7, and the experimental results of references 3 to 6 is not unex- 
pected since equation (3) is valid only in the temperature region where repulsion forces 
completely dominate. Also shown in figure 1 are values of B(T) obtained from the rela- 
tion of reference 2 

B(T) ~ (1.3436 X 10-2)(15.8922 - In T) 3 - 4.39 exp [^(2.4177 X 10-3) t] (6) 

This relation is a modification of equation (3); the addition of a negative exponential term 
makes it agree more closely with the experimental results of references 3 to 6 for tem- 
peratures down to 300° K. Following the example of reference 2, equation (3) was modi- 
fied in a similar manner to yield the relation 

B(T) * (1.3436 X 10- 2 )(15.8922 - In T) 3 - 8.04 exp -(3.7156 X 10" 3 )t] (7) 

Equation (7) extends the temperature range of equation (3) down to the region of maximum 
B(T), which occurs at about 200° K (see fig. 1) and will be used herein. The exponential 
term of equation (7) is negligible for temperatures greater than about 2000° K. The third 
and fourth virial coefficients agree reasonably well with experimental data for tempera- 
tures as low as 300° K (ref. 7). According to reference 7, the error incurred in using 
equation (3), and hence in using equation (7), may be as high as 5 to 10 percent for the 
present temperature range whereas errors incurred in using equations (4) and (5) may be 
as high as 25 percent. However, the virial coefficients B(T), C(T), and D(T) repre- 
sent perturbations in the ideal- gas behavior of helium to the first, second, and third order, 
respectively. Hence, as pointed out in reference 7, relatively large errors in the virial 
coefficients represent only small errors in the parameters of interest. 

Following the method of development of reference 1, the expression for enthalpy is 
found by considering enthalpy as a function of density and temperature 
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dh =(fW(w) p dT 


( 8 ) 


Introducing the combined first and second law equation of thermodynamics in the form 

dh = T ds + i dp (9) 

which can be rewritten for a constant temperature as 


and for a constant density as 




l/V 

P\ 0 P, 


T 





equation (8) may be expressed as 


dh = 


t/£p\ 

pA^Jp 



+ 




(10) 


(ID 


( 12 ) 


by noting, from the Maxwell relation, that 



W£p\ 

p 2 W p 


and 



(13) 


( 14 ) 


The pressure p is eliminated from equation (12) by use of equation (1). Direct differen- 
tiation of equation (1) gives 



pR < 1 + p 


B(T) + T 


dB(T) 

dT 



(15) 


and 


= RT 1 + 2pB(T) + 3p 2 C(T) + 4 p 
v 9 P/T L 



(16) 


These two equations are now substituted into equation (12), and equation (12), which is an 
exact differential, may be integrated over a two-step path by first holding p constant 
and integrating with respect to T and then holding T constant and integrating with 
respect to p. When integrating from a reference temperature to a desired temperature 
for a constant density, the density level chosen can be made as small as desired (p - 0) 
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3 

so that the value of c v in equation (12) may be set equal to its ideal-gas value of - R 

and f if— \ dT = RT. The second integration is performed for a constant desired 
J P \9T/p 

temperature from this very low reference density up to the desired density. The refer- 
ence value of the enthalpy (constant of integration) is eliminated by the reference- 
temperature term since the reference conditions correspond to ideal-gas behavior; that 
is, h = - RT for p — 0. Upon integration of equation (12), the expression for enthalpy 
as a function of p and T becomes 

h = RT<J| + p B(T) - T + Pi 2C(T) - T ^ 3D(T) - T ^2- J (17) 

where the virial coefficients B(T), C(T), and D(T) are given in equations (7), (4), 
and (5), respectively. Since the virial coefficients are functions only of the temperature, 
the derivatives as required in equation (17) can be obtained by direct differentiation. 

The expression for entropy as a function of p and T is found by considering 
which can be rewritten, by using equations (13) and (14), as 

ds = -^) p dp + T dT 


(18) 

(19) 


The pressure dependence is eliminated by using equation (15). By integrating over the 
same two paths for which the enthalpy was determined and noting that the reference value 
of the entropy (constant of integration) is eliminated at the reference conditions, the 
expression for entropy as a function of p and T becomes 



The expression for c v can be determined from equations (14) and (20) and is 

; T . f B -.WlfL jgill » T gjgail + gli jSM + T + fiffc jgg) + T fowl! 

v 2 | dT dT^ 2 dT dT^ 8 dT dT^ ( 

(21: 

The egression for the speed of sound is found by beginning with the relations 

a 2 = /M .!%(*) (22 

\®P/S c V\ 9 P/t 


fp _ _ jrM M . 1 

c v " c v p^ \ 9T /p\ aT /p 


(23) 
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Multiplying equation (23) by equation (24) and substituting the product into equation (22) 
yields / 


/ 9 P^ , T | 

^9PV 

o 2 

V9T ), 


(25) 


where the two quantities, 
respectively. 




, are given by equations (15) and (16), 
T 


The expression for c p can be found from equations (23) and (24) in the form 

,2 


c p - c v + 


n 2 \3T 


P 


9p\ 

9p/rp 


(26) 


where (3p/9T)p and (3p/ 9 P)x are given by equations (15) and (16), respectively, and 
c v has been determined previously from equation (21). 


CALCULATION OF FLOW PARAMETERS 


The previously derived equations for enthalpy, entropy, speed of sound, and specific 
heats are used, in conjunction with the equation of state, to provide expressions for deter- 
mining reservoir, nozzle-throat, test-section free-stream, and test-section normal-shock 
conditions in real helium for reservoir temperatures from 200° to 15 000° K and reservoir 
pressures from 1 to 3600 atmospheres (1 atmosphere = 0.101 MN/m 2 ). 

The following procedure for determining various thermodynamic and gas dynamic 
flow parameters was derived primarily to meet the needs of hotshot-type wind tunnels 
where, in many instances, the only direct measurements made to define the flow conditions 
are reservoir pressure after arc discharge and test- section stagnation pressure behind a 
normal shock. These measured pressures and the initial reservoir density before arc 
discharge are the basic input data. This procedure is based on the assumptions of the 
existence of uniform density throughout the reservoir during a tunnel test and an isen- 
tropic, one-dimensional flow expansion. The assumption of uniform reservoir density can 
be avoided by providing an additional measured quantity as input to the calculation proce- 
dure. Hence, three possible measured third inputs (free-stream density p^, free-stream 
velocity Uoo, and reservoir temperature T t ^ are also considered in the following cal- 
culation procedure. 



The following derived method for calculating reservoir, nozzle-throat, test- section 
free- stream, and test- section normal- shock conditions was programed for an electronic 
data processing system. The Newton- Raphson iteration method was employed to solve 
single equations and simultaneous equations where the initial conditions required for this 
iteration method were those for an ideal gas. 

Before the discussion of the calculation procedure for determining various flow 
parameters, a brief description of the Langley hotshot tunnel is given. The major com- 
ponents of this facility include a capacitor bank, an arc chamber incorporating coaxial 
electrodes, a conical nozzle and test section, and a vacuum reservoir. The operating 
procedure begins by charging (pressurizing) the arc chamber at ambient temperature with 
the test gas (helium or nitrogen). A quantity of stored electrical energy in the capacitor 
bank is then discharged across the coaxial electrodes. The resulting increase in pres- 
sure and temperature of the gas within the arc chamber results in the rupture of a dia- 
phragm between the arc chamber and evacuated nozzle. The gas then expands through the 
conical nozzle into the vacuum reservoir. A more detailed description of this facility is 
presented in reference 8. 

Initial Reservoir Conditions 

The initial reservoir density before arc discharge is found by using equation (1) in 
the form 

= p c + P c 2 B(T c ) + P c 3 C(T c ) + P c 4 D(T c ) (27) 

where the charge pressure p c is a measured quantity and the charge temperature T c 
has been experimentally found to become ambient shortly after charging the reservoir. 

The parameters B, C, and D are the virial coefficients (eqs. (7), (4), and (5)) and are 
functions only of charge temperature. 

Reservoir Conditions After Arc Discharge 

At a tunnel run time of zero (immediately after arc discharge), the reservoir density 
is equal to the initial reservoir charge density, that is, (p^ i) t _q = P c - K is assumed 
that no density gradients exist in the reservoir at any time after arc discharge. With the 
reservoir density and measured reservoir pressure known, the reservoir temperature 
for r = 0 can be calculated by using equation (1), in which the virial coefficients are 
functions of the reservoir temperature for r = 0. Having determined the reservoir den- 
sity and temperature, the reservoir compressibility factor, enthalpy, entropy, specific 
heat at constant volume, speed of sound, and specific heat at constant pressure can be cal- 
culated directly from equations (2), (17), (20), (21), (25), and (26), respectively. 

For r > 0, the time- monitored reservoir pressure is known, but the stagnation 
density has decreased from the value for t = 0 as a result of the mass flow from the 
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reservoir. To account for this time variation in stagnation density, it is necessary to cal- 
culate the mass flow from the reservoir (discussed in the following section) and adjust 
the reservoir density accordingly. 


Nozzle- Throat Conditions 


As mentioned previously, it is necessary to determine the mass flow from the 
reservoir through the nozzle throat as a function of tunnel run time. The relationships 
necessary for determining the mass flow from the reservoir through the nozzle throat 
are derived by assuming isentropic flow 

s*(p*,T*) = s tjl (28) 

and noting that the Mach number at the nozzle throat M* is unity; that is, the speed of 
sound at the nozzle throat is equal to the velocity at the nozzle throat (a* = U*). The 
energy equation 

U 2 = 2(h t>1 - h) (29) 

is considered in conjunction with a* = U* to yield the expression 

h*(p*,T*) + i[a*(p*,T*)] 2 = h t> ! (30) 

where the reservoir enthalpy h+ i has been calculated previously, h* is given by equa- 
tion (17) as a function of p* and T*, and ( a *)* is given by equation (25) also as a 
function of p* and T*. Since s t j has been calculated previously and s* is given 
by equation (20) as a function of p* and T*, equations (28) and (30) can be solved simul- 
taneously for p* and T*. With p* and T* determined, h* may be calculated 
directly from equation (17) and U* from equation (29). 

The reservoir density can now be adjusted (adjustment is required due to mass flow 
from the reservoir with tunnel run time) by the equation 


( p m) t "KOt 


present 


previous 


0 *TJ* 7T(dTAT 

H 4V 


where the nozzle-throat diameter d* and reservoir volume Y are geometric input 
parameters and Ar represents the desired time increment. Depending on the rate of 
mass flow from the reservoir through the nozzle throat, the time increments of equa- 
tion (31) should be chosen small enough so that a further decrease in At does not result 
in significant variations in p t j for all r (larger mass flow rates require smaller 
time increments). 

The determination of the reservoir density, hence reservoir temperature, enthalpy, 
and so forth, as a function of tunnel run time is best illustrated by an example. Assume 
that p t j has been measured at 0.005- sec intervals, beginning immediately after arc 
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discharge (r = 0), so that rg = 0, t\ = 0.005 sec, r 2 = 0.010 sec, .... For tq = 0, 

Pt 1 = p c , thereby permitting calculation of h^j, s^j_, and so forth, at tq. 

Since only the measured quantity p^. j is known at tj, p^j at tj is found by sub- 
tracting the mass flow from the reservoir over the time increment t\ - Tq from the 
previously calculated value of p ^ j at tq. This is accomplished by first calculating 
p* and T* from the simultaneous solution of equations (28) and (30) where ht ? i and 
s^ correspond to tq. The nozzle- throat velocity U* is calculated from equation (29), 
where h^ and h* correspond to tq, and equation (31) becomes 


Kl) T f ■Ki)t 0 - (p * D * ) t 0 ;1:L ^ 


ir(d*) 2 (r x - t 0 ) 


This value of (p^ is the first approximation (indicated by superscript (1)) of p^ ^ 

at t-^. By using ^p t i) T ^^ in conjunction with the measured value of (Pt 1 ) , the 
parameters T^ j, hj. j., Sj. j, p*, and U* can be calculated at and a new value 

\2 


of p ^ at tj can be found from 




( 1 ) 


7r(d*) 


( T 1 ~ T o) 


where 


(Pt, l) T ^ is the second approximation of P t ^ at tj and (p*U*) T ^ 


( 1 ) 


corre- 


sponds to the first approximation of p t j at tj obtained from the preceding equation. 


If this new value of 


Pt, 1 at rj ^that is, ^p^ ^ ( 2 )j is within the desired accuracy 


(0.01 percent for the present study) of the previous value of P t ^ at tj (that is, 

(1)\ ’ ' 

^Pt, i) T ^ J the calculated parameters corresponding to the new value of p ^ j at Tj are 


considered satisfactory; if not, this procedure must be repeated until the desired accuracy 
is obtained. The nozzle-throat pressure p* at tj can be calculated directly from 
equation (1) since p* and T* at tj are now known. 


Free- Stream Conditions 


For Moo > 10, the free-stream enthalpy is less than about 3 percent of the reser- 
voir enthalpy. Therefore, the free-stream velocity can be approximated by allowing 
hoo = 0 in equation (29), which gives 

U 0O 2 = 2h t;1 (32) 

The ratio of the free-stream dynamic pressure to reservoir pressure for isentropic 
flow of an ideal gas is given in reference 9 as 
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(33) 


^= I m “ 2 ( 1 + V m " 2 ) ’'- 1 

Also from reference 9, the ratio of total pressure behind a normal shock to reser- 
voir pressure for adiabatic flow of an ideal gas is given by 


P t ,2 

(y + DM, 2 

y-l 

y+1 

p t,r 

(y - 1 )mJ + 2_ 


2yMj - (y - 1) 


1 

y-l 


(34) 


1 2 

Dividing equation (33) by equation (34) and noting that q^ = results in the fol- 

lowing equation for free-stream density: 


p JhAH 

oo 


00 u, 2 


r 


-l 


y 


y+ 1 

r-i/ 


y+1 


1 

k M 2 


1 

y-l 


(35) 


Since M^ 2 » 1 for the present study, equation (35) can be simplified to 

p„ “ r(,A| 

Uco 2 


(36) 


where 


r(y) = 2 


1 y+1 



and r(y) is equal to 1.1347 for ideal helium (y = 1.6667). A superscript (~) on p ro 
indicates an approximate value since equation (36) was obtained from ideal-gas relations. 
Equation (36) thus provides an approximation of the free-stream density based on the mea- 
sured total pressure behind a normal shock and estimated free-stream velocity as deter- 
mined from equation (32). 

From the assumption s t j = Soo, and with known, can be calculated 

numerically from equation (20)’. Free-stream pressure p^ and free-stream enthalpy 
hoc can then be determined directly from equations (1) and (17), respectively. 

Since U^ represents an approximation of the free-stream velocity due to the 
omission of the free-stream enthalpy contribution, free-stream velocity can now be cal- 
culated from equation (29) in the form 

U„ 2 = 2(h M - 

By using this value of Uoo, a- new value of p, Tj and hence of Too and can be cal- 
culated, thus leading to a second new value of Uoo- This procedure is repeated until the 
desired accuracy between successive values of Uoo (0.01 percent for the present study) 
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is obtained. The free- stream parameters (p^Too) corresponding to the final value of 
free-stream velocity in this iteration can now be used to calculate the free-stream speed 
of sound aoo from equation (25). The free-stream Mach number Moo is simply found 
by dividing Uoo by aoo, and the free-stream dynamic pressure is found from 

1~ p 

Q .00 “ 2 /3 °o^ 00 • The values of Z^, c Vj00 , and Cp >00 are calculated directly from equa- 
tions (2), (21), and (26), respectively, and the free-stream ratio of specific heats is found 
c 

from y = The free-stream coefficient of viscosity /ioo is given to within 1 per- 

cent for the temperature range 4° to 1100° K by the empirical expression (ref. 10) 


Moo = (1.2548 X 10-6) Too 0 - 647 (37) 

where ^ is in moles/cm-sec for Too in °K. The experimental viscosity results 
collected in reference 10 show that for a temperature of 1.64° K, the value of q calcu- 
lated from equation (37) is 1.255 times the experimental value. However, it should be 
pointed out that this is the only available experimental result for temperatures less than 
4° K. In the absence of stronger evidence of error, equation (37) is employed in this 
study for Too < 4° K. The free-stream Reynolds number, per unit length, is found by 
dividing the product p M Uoo by ju^. The ratio of effective nozzle-throat area to effec- 
tive nozzle-test- section area is obtained from the equation of continuity 

A* _ Poo^oo 

“a " *77*" 

A p U 


(38) 


In order to calculate these free-stream conditions, the free-stream velocity and 
density were approximated from equations (32) and (36), respectively. Although the 
inaccuracy in approximating free-stream velocity was eliminated in this section by an 
iterative procedure, the free-stream density is still an approximation since equation (36) 
was derived by using ideal-gas relations. The validity of this value of is examined 
in a subsequent section. 

In the calculation of these free-stream conditions, the real-helium expressions 
developed in a previous section "Basic Thermodynamic Relations” were used. For the 
present range of reservoir conditions and Mach number, the free-stream temperature 
range is approximately 1.0° to 440° K and the maximum free-stream density is approxi- 
mately 1.2 x 10“ 4 mole/cm^. Due to this low range of free-stream density, terms 
appearing in equation (2) which involve C(T) and D(T) are essentially negligible in the 
free stream. The use of equation (7) for determining B(T) at free-stream temperatures 
less than 200° K is questionable. (See fig. 1.) For the most unfavorable conditions, the 
use of equation (7) results in a value of Z M which is approximately 1.003 times the more 
correct value obtained from the experimental values of B(T) tabulated in reference 10 
for very low temperatures. By allowing B(Too) =0 for Too < 200° K, this discrepancy 
in Zqq is reduced to approximately 1.001; therefore, in the present method, B(Too), 
C(Too), and D(Too) were set equal to zero when Too < 200° K. 
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An item of concern in the free- stream flow of a gas is that of flow condensation. 

To determine whether helium condensation might be present in the free stream, the fol- 
lowing empirical expressions from reference 10 are included in the present method. For 
Too £ 2.190° K 

log 10 Peon = " ^7 + 2 - 208 log 10 To ° ' °- 664 (39a) 

and for Too < 2.190° K 


logjo P con = - ^ + 0.922 log 10 T„ + 0.154 (39b) 

where p CQn is in atmospheres for Too in °K. If the free-stream pressure p^ is 
less than the calculated condensation pressure P con , condensation will not be expected 
to occur; however, if p^ ^ p CQn , condensation may occur. 


Static Conditions Behind Normal Shock 


The conservation relations for mass, momentum, and energy across a normal shock 

are 


PooU°° - P2 U 2 

(40) 

oo U oo 2 = P 2 + p 2 U 2 2 

(41) 


and 

= ho. + \ Uoo 2 = h 2 + \ U 2 2 


(42) 


From equations (1), (40), and (41) the expression 

P. + P„U„ 2 = P 2 RT 2 fl + P 2 b(t 2 ) + P 2 2 c(t 2 ) + P 2 3 d(t 2 ) 


(p»p»y 

+ p 2 


can be obtained. Similarly, from equations (17), (40), and (42) the expression 


bt,l ~ RT 2 < W + P 2 


p/ T \ T dB ( T 2) 

B V T 2) ' T 2 


2C(T 2 ) - 


dC 


2 dT 


^ 2 ) 


(43) 



(44) 


can be obtained. Equations (43) and (44) can be solved simultaneously for p 2 and Tg. 
With p 2 and T 2 known, the parameters p 2 , h 2 , s 2 , a 2 , c^ 2 , c p>2 , y 2 , and Z 2 
can be calculated directly. 
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At this point in the calculation procedure, the value of free- stream density used in 
the calculation of free-stream conditions and static conditions behind a normal shock is an 
approximate value since equation (36) was derived by using ideal-gas relations. It is now 
possible to examine the validity of this value of free-stream density. The variation in 
conditions from immediately downstream of the shock to the stagnation point is relatively 
small, and the level of density in this region is also relatively small for the present range 
of reservoir conditions and free-stream Mach number (approximately 4 times the maxi- 
mum free-stream density). Thus, this region may be considered to behave as an ideal 
gas and must satisfy the condition s 2 = s t 2 - Now p 2 can be calculated from the ideal- 
gas isentropic relation (ref. 9) 


- 1 


p 2 = p t,2 t 1 + yjL r~ M 2 2 ; 


^2 


-1 


(45) 


where p^. 2 is the measured stagnation pressure behind a normal shock. This value of 
p 2 is compared with the value of p 2 obtained previously, and if these values are not 
within the desired accuracy (0.01 percent for the present study), a new value of is 
determined from the relation 

= (46 > 

where p 2 ’ denotes the value of p 2 obtained from equation (45). The calculation pro- 
cedure, beginning immediately after equation (36), where = (Poo)new’ * s re P ea ^ e< ^ until 
the desired accuracy is obtained. 


Stagnation Conditions Behind Normal Shock 


The stagnation conditions behind a normal shock must satisfy the condition 
2 = h t -p Thus, equation (17) can be expressed as a function of 2 and T t 2 . Since 
Pj. 2 is a measured quantity, equation (1) can also be expressed as a function of 2 and 
Tt 2- After p t 2 and T t 2 have been obtained from the simultaneous solution of equa- 
tions (1) and (17’), the parameters s t 2 , (c v ) t 2 , 2 , y t 2 , and Z^ 2 can be calcu- 

lated directly. Equation (37) can be used to calculate p T 2 for values of T t 2 less 
than 1100° K. To obtain an expression for 2 for values of T^ 2 greater than 
1100° K, a manual adjustment was made on equation (37) to obtain better agreement 
with results of reference 7; the resulting equation is 


M 


t,2 


(l.2548 X 10" 6 )T t>2 - 647 exp (6.04 X lO' 5 )^ 2 - (o.23 X 10“ 8 )T t 


(47) 
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where 2 is in moles/cm-sec for in °K. For values of T^2 less than 

1625° K, equation (37) predicts ju. t 2 within 4 percent of the results of reference 7, 
whereas for values of T t 2 between 1625° and 11 000° K, equation (47) yields values of 
2 within 4 percent. 


Procedure Using an Additional Experimental Input 

The previously discussed procedure for determining various flow parameters was 
derived for the case when only reservoir pressure and stagnation pressure behind a nor- 
mal shock are measured to define the flow. This procedure was based on the assumption 
that uniform density existed throughout the reservoir. The validity of this assumption is 
questionable in explosive arc-heated facilities, and this assumption can be avoided by mea- 
suring one other flow parameter in conjunction with the measured pressures. (It is 
assumed that simultaneous measurements of reservoir pressure, stagnation pressure 
behind a normal shock, and the third flow parameter are made throughout the test.) 

Free- stream density .- Successful free- stream density measurements have recently 
been made in the Langley hotshot tunnel by using an electron beam in nitrogen flow. If 
the electron-beam technique proves to be applicable for measuring free- stream densities 
in helium, this measurement can be utilized as an additional input. The calculation pro- 
cedure for determining the various flow parameters of interest when free- stream density 
is a measured quantity begins by combining equations (32) and (36) to form 


V 




(48) 


where p^. 2 is also a measured quantity and hoo has been neglected. This calculated 
value of reservoir enthalpy can be used in equation (17), thereby expressing equation (17) 
in terms of reservoir density and temperature. Since reservoir pressure is a measured 
quantity, equation (1) can also be expressed in terms of reservoir density and tempera- 
ture. Equations (1) and (17) can then be solved simultaneously for p ^ j and T^j.. When 
these reservoir parameters and, of course, p t ^ are known, the remaining phase of the 
calculation procedure is the same as discussed previously, except that the mass-flow 
iterative procedure (eq. (31)) is no longer necessary. Since equation (48) is only an 
approximation of the reservoir stagnation enthalpy, the final calculated value of p x is 
compared with the measured value of p^. If these values are not within the desired 
accuracy (0.01 percent for the present study), the reservoir density corresponding to the 
final calculated free- stream density is adjusted by use of the relation 



( 49 ) 
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This new value of reservoir density is introduced into equation (1) in order to calculate 
a new value of reservoir temperature. The calculation procedure is repeated with these 
new values of p t ^ and T t and a second comparison of the final calculated free- 
stream density, corresponding to these new values of p t ^ and T t is made with the 
measured free- stream density. This procedure is repeated until the desired accuracy 
is obtained. 

Free-stream velocity.- Measurements of free-stream velocity have been performed 
in the Langley hotshot tunnel by photographically observing the travel of spark discharges 
in the nitrogen flow stream. It is believed that this technique may also be applicable for 
determining free-stream velocities in helium. The calculation procedure for determining 
the various flow parameters of interest when free-stream velocity is a measured quantity 
begins by estimating the reservoir enthalpy from equation (32) where hoo is neglected 
for the time being. With h t 1 (approximately) and p t x known, equations (1) and (17) 
can be solved simultaneously for p t 1 and T^i- The’remaining reservoir and nozzle- 
throat conditions can be calculated as discussed previously, again noting that the mass- 
flow iterative procedure (eq. (31)) is no longer necessary. The final value of free-stream 
velocity calculated from equation (29), where the contribution of hoo is included, is com- 
pared with the measured value of free-stream velocity. If these values are not within 
the desired accuracy (0.01 percent for the present study), the left side of equation (17) is 
(tt V ^ 

set equal to > + h M , where the value of h^ is that corresponding to the final calcu- 

lated free-stream velocity. Thus, a second final value of free-stream velocity can be cal- 
culated and compared with the measured free-stream velocity. This procedure is 
repeated until the desired accuracy is obtained. 

Reservoir temperature .- The present technique for measuring reservoir tempera- 
ture in hotshot tunnels is not fully developed because of the severity of the arc-chamber 
conditions. Nevertheless, the continuous measurement of reservoir temperature is also 
considered as an additional input. The calculation procedure is then very similar to that 
discussed previously when only p t j and p t 2 are measured inputs, the only difference 
being that the reservoir density can now be calculated directly by the use of equation (1). 
Therefore, the mass-flow iterative procedure (eq. (31)) is no longer necessary. 

Hypothetical case using free-stream Mach number as an input. - Calculations of a 
number of hypersonic nozzle flow parameters for real helium have been performed and 
presented in reference 1 in the form of simple correction factors for reservoir pres- 
sures to 3400 atmospheres, temperatures to 5800° K, and free-stream Mach numbers 
to 100. The results of reference 1 were obtained by using the Beattie-Bridgeman equa- 
tion of state, for which the constants used were determined from experimental data for 
temperatures to only 590° K. In order to obtain more accurate real-helium effects at 
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temperatures above 590° K, the present method based on the virial equation of state was 
used to determine the flow properties of real helium by using free- stream Mach number 
as an input along with reservoir pressure and temperature. The purpose of this hypo- 
thetical calculation procedure is to obtain simple correction factors in the form of ratios 
of a real-helium value of a particular flow parameter to its ideal-helium value and to 
present these correction factors in working charts. 


The calculation procedure for determining reservoir and nozzle-throat conditions 
when reservoir pressure and temperature are inputs has been discussed previously. The 
free-stream conditions are determined by first estimating Uoo by using equation (32), 
where h^ is neglected, and considering the relation 


2 U„ 


a oo(^oo> r ^ 0 °) _ 




(50) 


where a^ is given by equation (25) and the right side of equation (50) is known. From 
the assumption of isentropic flow (s o0 (p 0O ,'To^ = s t equations (20) and (50) can be solved 
simultaneously for and T M ; thereby, the use of equation (36) and hence of equa- 
tions (45) and (46) is avoided. The free-stream enthalpy can be calculated directly from 
equation (17) and a new value of Uoo can be found from equation (29). This new value of 
Uoo is compared with the previously calculated value of U^. If these values of U^ are 
not within the desired accuracy of each other (0.01 percent for the present study), the new 
value of Uoo is placed into equation (50) and this procedure is repeated until the desired 
accuracy between successive values of U^ is obtained. With and Too known, 

the various free-stream parameters of interest can be calculated. The static conditions 
behind a normal shock can be determined as discussed previously. The stagnation condi- 
tions behind a normal shock must satisfy the conditions h^. 2 (p^ 2’^t 2 ) = ^t 1 anc * 
s t 2( i °t 2’^t 2 ) = s 2’ w ^ ere 1 an< * s 2 h ave t> een determined previously. Thus, equa- 
tions (17) and (20) can be solved simultaneously for p^ 2 and T^ 2 and p t 2 can be 
determined from equation (1). The results of this hypothetical calculation procedure are 
restricted to reservoir densities less than 0.0173 mole/cnr*, the critical density of helium, 
and to uncondensed free-stream flow as defined by equation (39). 


RESULTS AND DISCUSSION 

The results of the calculations for helium based on the virial form of the equation 
of state where free-stream Mach number was used as an input along with reservoir pres- 
sure and temperature are shown in figures 2, 3, 4, and 5. Various reservoir and nozzle- 
throat thermodynamic parameters are plotted in figures 2 and 3, respectively, as a func- 
tion of reservoir pressure to 3600 atmospheres for reservoir temperatures ranging from 
200° to 15 000° K. The ideal-gas value of a particular reservoir or nozzle-throat 
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thermodynamic parameter is shown in figures 2 and 3 where appropriate. (The ideal-gas 
expressions of ref. 11 are shown in fig. 3.) The results of figures 2 and 3 are restricted 
to reservoir densities less than the critical density of helium. 


Figure 4 shows the results of the calculations in the form of correction factors to 
be used in conjunction with the ideal-gas parameters for helium presented in refer- 
ences 12 and 13. These correction factors are ratios of real-to- ideal flow parameters 
and are plotted in figure 4 as functions of reservoir pressure to 3600 atmospheres for 
reservoir temperatures ranging from 200° to 15 000° K and a range of Mach numbers 
from 10 to 100. The Mach number ratio of figure 4(a) represents the real-helium, or 
input, Mach number divided by the ideal- helium Mach number as determined from equa- 
tion (34), where p^. 2fit \ is a real-helium quantity ^ ^ is an input parameter and 
Pj. g corresponds to input Mach number J and y - 5/3. As in reference 1, the real- 
helium flow parameters of figures 4(b) to 4(f) are divided by the corresponding ideal- gas 
flow parameter (obtained from refs. 12 or 13) for the same value of Moo to give the cor- 
rection factors. The results of figure 4 are restricted to reservoir densities less than 
the critical density of helium and to Mach numbers corresponding to uncondensed flow. 

The use of the correction factors presented in figure 4 is essentially the same as 
discussed in reference 1 except that no iterative procedure is required for the determina- 
tion of the real- helium Mach number. For example, assume that a hypersonic helium 

tunnel operates with p. = 1000 atm, T t j = 900° K, and Pj 2 ~ °-9533 atm. From 

, v ’ ’ M_ ’ P t o 

figure 4(a), these conditions correspond to ^ y = 1.0429, and for = 0.9533 x 10"°, 

reference 12 yields (M 00 ) i = 28.77; therefore, M^ is determined by multiplying the 
ideal-gas value of Mach number as found in reference 12 by the correction factor 1.0429 
and is found to be equal to 30. Once the value of M^ has been determined, other cor- 
responding ideal- gas flow parameters can be found from the tables of reference 12 
(1 = Moo = 40) or reference 13 (40 § Moo = 100). For instance, for Moo = 30, reference 12 

shows that ( 5 —'] = 0.6362 x 10 "® and figure 4(b) shows that for P^. j = 1000 atm, 

V t, 1/ x ’ 

T t i = 900° K, and M^ = 30, ^^^-= 1.1330; Poo fit 1 is found, by multiplying the 


1-6. 


p _ = 0.7208 x 10 -3 atm. 


ideal-gas value by the correction factor, to be 0.7208 X 10 , Coo 

The other flow parameters are determined in the same manner. For Mach numbers 
between 10 and 30, a Mach number dependency is observed for most of the correction 
factors; however, for Mach numbers greater than 30, the correction factors are essen- 
tially independent of Mach number. The correction factors for P 2 /P 00 ’ ^2 ;/°oo» ^ 2 /^°°’ 
p t 2 /P 00 ’ an( * m 2 are n °t Presented since they were all within 0.3 percent of unity for 
a Mach number of 10 and within 0.05 percent of unity for a Mach number of 20 or greater 
for the entire range of reservoir pressure and temperature. The correction factors for 
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qco /p t 1 and p t 2 yp t 1 were within 0.05 percent of the correction factor for p^/Pj. j 

for given p t ^ T t and M m for the entire range of calculation and hence were con- 
sidered to be’ equal (fig. 4(b)). To avoid the cross-plotting generally associated with 
working charts, the correction factors of figures 4(a) to 4(d) were cast into the form of 
simple equations which are presented in the appendix. 

In figure 5, the compressibility factor Z t ± is plotted as a function of reservoir 
pressure for various reservoir temperatures. The solid lines denote the results of the 
present calculations based on the virial equation of state, whereas the dash lines denote 
the results of the calculation procedure of reference 1 which was based on the Beattie- 
Bridgeman equation of state. The compressibility factor was selected to make a direct 
comparison between the present calculation procedure and that of reference 1 since it is 
the governing parameter for real-gas effects in helium. Figure 5 shows that the com- 
pressibility factors calculated in reference 1 were always greater than those of the pres- 
ent calculations for a given pressure and temperature. Because of the larger calculated 
compressibility factors of reference 1 (which imply a greater departure from ideal- gas 
behavior) the correction factors of reference 1 tend to overcorrect the ideal-gas param- 
eters in comparison with the present correction factors. For example, for 
p t j = 2000 atm and T t j = 1367° K, a 5.1 percent discrepancy in Z t> j corresponds to 
the correction factors for p^/Pj. v P o0 /P i v and Too/T t>1 of reference 1 being approxi- 
mately 7.3 percent, 11.4 percent, ’and 1.9 percent greater, respectively, than those of the 
present study. The existing discrepancies between the results of the present study and 
those of reference 1 are not unexpected since, as pointed out in reference 1, the results 
of reference 1 were obtained by using constants for the Beattie-Bridgeman equation of 
state which were determined from experimental data for temperatures to only 590° K and 
thus should be treated as approximations for T^ ^ > 590° K. Hence, it is believed that a 
more accurate description of real-helium effects is obtained from the present calculation 
procedure based on the virial equation of state. 

As discussed previously, the expression for free- stream density (eq. (36)) was 
obtained from ideal-gas relations, thereby making it necessary to examine its validity. 

This was performed by assuming that the region from immediately downstream of the 
shock to the stagnation point behaves as an ideal gas and thus allows equation (45) to be 
employed. The results of the calculations of the hypothetical case (using free- stream 
Mach number as input), which avoid the use of equation (36), showed that p 2 determined 
from equation (45) (denoted p 2 ') was within 0.04 percent of p 2 determined from equa- 
tion (1) at Moo = 10 ^ within 0.005 percent at Moo ^ 20 for the present range of reser- 
voir pressure and temperature. This supports the validity of equation (45) and hence the 
use of equation (46) to determine the free- stream density. Calculations for a number of 
tunnel tests at Mach numbers from approximately 30 to 80 showed that equation (36) is 
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very accurate, yielding results within 0.3 percent of the final value of obtained from 
equation (46) at ~ 30, with the accuracy improving with increasing Mach number. 

The accuracy of equation (36) can be attributed to the fact that the ideal-gas relations 
(eqs. (33) and (34)) used to derive equation (36) have real-helium correction factors 
within 0.05 percent of one another, for given p^. j, Tf. j, and M m , over the present 
range of calculation. 

CONCLUDING REMARKS 

The virial form of the equation of state was used to derive a method for deter- 
mining reservoir, nozzle-throat, free- stream, and normal- shock conditions for real- 
helium flow. This method was derived primarily to meet the needs of hotshot-type tun- 
nels, where the basic input data are measured reservoir pressure and stagnation pressure 
behind a normal shock, and was based on the assumption of uniform density throughout the 
reservoir during a tunnel test. The method was derived with the capability of using an 
additional measurement of free- stream density, free- stream velocity, or reservoir tem- 
perature; this additional measurement as input allows the assumption of uniform density 
to be avoided. The method was applied for reservoir pressures to 3600 atmospheres, 
temperatures to 15 000° K, and a range of free-stream Mach numbers from 10 to 100, 
and the results are presented in the form of nondimensionalized reservoir thermodynamic 
and nozzle-throat parameters and as simple correction factors which may be applied to 
ideal-helium flow parameters. The present correction factors were compared with the 
correction factors of NASA TN D-1632, which were obtained by using the Beattie- 
Bridgeman equation of state where the constants used were determined from experimental 
data for temperatures to only 590° K. This comparison showed that the correction fac- 
tors of NASA TN D-1632 tend to overcorrect the ideal-helium flow parameters. The 
present calculations are restricted to reservoir densities less than the critical density of 
helium and to theoretically uncondensed free-stream flow. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 24, 1968, 

124-07-02-54-23. 
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APPENDIX 


SIMPLIFIED METHOD FOR DETERMINING REAL- HELIUM 
HYPERSONIC FLOW PARAMETERS 


To provide a rapid determination of various real-helium flow parameters, the cor- 
rection factors of figures 4(a) to 4(d) were cast into the form of simple equations, and 
thereby the cross-plotting generally associated with working charts such as figure 4 was 
avoided. 


Except for Mach numbers between 10 and 30, the correction factors were essen- 
tially independent of Mach number. Consequently, equations were obtained for each cor- 
rection factor in terms of reservoir temperature and pressure for Moo = 30 and 
300° K = T^ i = 15 000° K. The equations are 

(51) 


(52) 


(53) 


(54) 

where p t j is in atmospheres and T t j is in °K. Equations (52) to (54) may be applied 
for Mach numbers between 20 and 30 with a small uncertainty (less than 0.2 percent at 
M^ = 20). Use of these equations for Mach numbers less than 20 will result in correc- 
tion factors that are as much as 1 percent too large for equations (52) and (53) and 
0.5 percent for equation (54) for the present range of reservoir conditions. 
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Figure 1.- Second viriai coefficient for helium as a function of temperature. 
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Figure 2.- Reservoir thermodynamic parameters for helium as a function of reservoir pressure for various reservoir temperatures. 
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Figure 5.- Comparison of compressibility factors obtained with present calculation procedure and those of reference 1. 



National Aeronautics and Space Administration 
Washington, D. C. 20546 

official business FIRST CLASS MAIL 


POSTAGE and fees paid 
NATIONAL AERONAUTICS AND i 
SPACE ADMINISTRATION f 


t 


POSTMASTER: 


If Undeliverable ( Section 158 
Postal Manual ) Do Not Return J 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute . . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 

NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 

TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES : Information less broad 
in scope but nevertheless of importance as a 
. contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 

Details on the availability of these publications may be obtained from: 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



